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ABSTRACT 

We present the results of a search for galaxy alignments in 12 galaxy clusters at 
z > 0.5, a statistically complete subset of the very X-ray luminous clusters from the 
MAssive Cluster Survey (MACS). Using high-quality images taken with the Hubble 
Space Telescope (HST) that render measurement errors negligible, we hnd no radial 
galaxy alignments within 500 kpc of the cluster centres for a sample of 545 spec- 
troscopically confirmed cluster members. A mild, but statistically insignificant trend 
favouring radial alignments is observed within a radius of 200 kpc and traced to galax- 
ies on the cluster red sequence. Our results for massive clusters at z > 0.5 are in stark 
contrast to the findings of previous studies which find highly significant radial align- 
ments of galaxies in nearby clusters at z ~ 0.1 out to at least half the virial radius 
using imaging data from the SDSS. The discrepancy becomes even more startling if 
radial alignment becomes more prevalent at decreasing clustercentric distance, as sug- 
gested by both our and previous work. We investigate and discuss potential causes for 
the disparity between our findings based on HST images of clusters at z > 0.5 and 
those obtained using groundbased images of systems at z ~ 0.1. We conclude that 
the most likely explanation is either dramatic evolution with redshift (in the sense 
that radial alignments are less pronounced in dynamically younger systems) or the 
presence of systematic biases in the analysis of SDSS imaging data that cause at least 
partly spurious alignment signals. 
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1 INTRODUCTION 

Galaxy and cluster orientations provide crucial insights into 
the processes governing cluster formation and evolution. 
On large scales, the distributions of galaxies in clusters are 
found to align with the directi on toward thei r neighbours 
within distances of tens of Mpc ((BiriEEcli 198^), a tendency 
that is explained by a preference for clust ering and mergin g 
along large-scale filamentary structures l|West et al.lll995T ). 
On smaller scales, intrinsic galaxy alignments generated dur- 
ing galaxy forma tion are observed on scales of a few Mpc 
i|Lee fc Penll2007l ') with different origins, however, for ellipti- 
cal and spiral galaxies. Similar to the physical mechanisms 
behind cluster alignments, anisotropic accretion of material 
onto galaxy halos along filamentar y structure s can result in 
the alignment of elliptical galaxies (West 1994). By contrast, 
the alignment of spiral galaxies is believed to be caused by 
correlations between t heir spin axes or iginating from initial 
large-scale tidal fields l|Pen et al.ll2000t ). 
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In dense cluster environments, extensive evidence is 
found of alignments between the brightes t cluster galax- 
ies (BCGs) and their host-cluster halos ([Binggelil 1 19821 ). 
an effect that is again caused by dynamical connections 
to large-scale filamentary structures. For cluster members 
other than the BCGs, any intrinsic alignment might be ex- 
pected to be destroyed by their gravitational interaction 
with the host-clu ster halos on a time scale of a few or- 
bits (|Couttslll996l ). However, using Sloan Digital Sky Sur- 
vey (SDSS) data in a study of a sample of 85 X- r ay se - 
lected clusters at 0.02 < z < 0.23, iPereira fc Kuhnl (|2005h 
(hereafter PK05) find a tendency for the major axes of clus- 
ter galaxies to be aligned with the radius vector pointing 
toward the cluster centre. They propose that this radial 
alignment is the result of tidal torquing exerted by the host 
clust e r halo (IPereira. Bryan, fc Gilll 120081 : IPereira fc Brvanl 
l2010h . lFaltenbacher et alj (|2007l ) (hereafter F07) confirm the 
radial alignment signal for a sample of galaxy groups at 
0.02 < z < 0.2 selected from the SDSS group catalogue, 
and further find a trend for the signal to be more significant 
near the cluster centre and for redder galaxies. 
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So far almost all searches for galaxy alignments have 
been conducted for relatively nearby clusters owing to the 
limited depth of readily available ground-based imaging 
data. If a study akin to PK05 could be conducted at sig- 
nificantly higher redshift, the presence or absence of galaxy 
alignments in distant clusters could shed light on the tem- 
poral evolution of the effect and thus on the dynamical 
timescales of the physical mechanisms at work. Unfortu- 
nately, measuring galaxy orientations from SDSS data is 
extremely challenging already at low to moderate redshifts 
due to the poor seeing and insufficient pixel sampling. High- 
quality images are thus essential in obtaining precise galaxy 
orientation measurements, specifically at high redshift. 

In this paper, we discuss our study of galaxy alignments 
in 12 X-ray luminous cl usters at z > 0.5 fro m the Massive 
Cluster Survey (MACS; lEbeling et all 120071 ') as a first step 
toward investigating the evolution of galaxy alignments with 
lookback time. All 12 clusters were observed with the Ad- 
vanced Camera for Surveys (ACS) aboard the Hubble Space 
Telescope (HST). We describe our target clusters, the HST 
observations, and image reduction procedures in §2. Clus- 
ter membership determination and galaxy shape measure- 
ments are described in §3. Results on galaxy alignments are 
presented in §4.1 and §4.2, and possible systematics are ex- 
amined in §4.3. We interpret our results and compare them 
with those of studies in nearby clusters in §5; a conclusion 
is given in §6. 

AB magnitudes are used throughout. We adopt the con- 
cordance ACDM cosmology with Ho = 70 km s~ Mpc~ , 
M = 0.3 and fl A = 0.7. 



2 TWELVE MACS CLUSTERS AT Z > 0.5 

MACS provides a statistically comple te sample of the mos t 
X-ray lumino u s clus ters at z > 0.3 (jEbeling et al.ll200l| y 
lEbeling et ail (|2007h present a subsample of the 12 most 
distant MACS clusters at z > 0.5, with 11 of them at 0.5 < 
z < 0.6 and MACS J0744.8+3927 at z = 0.698. Like all 
MACS cl usters, these system s are very massive (M ~ 10 14 — 
10 15 M Q ; iMantz et all bold ). Being purely X-ray selected, 
this sample is highly diverse with respect to optical ap- 
pearance and dynamical state. For example, MACS J0025.4- 
1222 consists of two substructures of nearly equal mass with 
the hot gas residing in between them, showing it to be a 
post- collision merger ( fBradac et al.l 120081 : iMann fc Ebelind 
|2011| ). Similarly, in-depth analyses of MACS J0717.5+3745 
and MACS J11 49.5+222 show highly complex s tructures 
and dynamics (|Ma et al.l [20091 ; ISmith et aLN2009h . On the 
other hand, MACS J1423. 8+2404 features a prominent cool 
core and a single domi nant cD galaxy, impl ying that it is 
close to fully virialized dLimousin et al.ll2010t ). 

2.1 HST Imaging 

All 12 target clusters were imaged with the ACS (Wide 
Field Channel) on HST (PID: 9722, 10703, 11560, PI: 
H. Ebeling; PID: 10493, PI: A. Gal- Yam). The field of 
view of ACS/WFC (3.5 x 3.5 arcmin 2 ) corresponds to ap- 
proximately 1.23 x 1.23 Mpc 2 at z~0.55 and is thus well 
matched to the physical size of the cluster cores. Each clus- 
ter field was imaged in the F555W and F814W filters for 



Table 1. Numbers of cluster members in the 12 MACS clusters 
within the ACS field of view for various subsamples: spectroscopi- 
cally confirmed (A), red-sequence selected (B), both spectroscop- 
ically and colour-selected (A+B) and spectroscopically but not 
colour selected (A-B). 



MACS name 




n(A) 


n(B) 


n(A4-B) 


n(A-B) 


MACSJ0018.54 


-1626 


100 


233 


62 


38 


MACSJ0025.4- 


-1222 


74 


223 


62 


12 


MACSJ0257.1- 


-2325 


50 


142 


33 


17 


MACSJ0454.1- 


-0300 


75 


159 


52 


23 


MACSJ0647.74 


-7015 


30 


71 


12 


18 


MACSJ0717.54 


-3745 


120 


316 


97 


23 


MACSJ0744.84 


-3927 


65 


160 


46 


19 


MACSJ0911.24 


-1746 


63 


162 


46 


17 


MACSJ1149.54 


-2223 


80 


267 


58 


22 


MACSJ1423.84 


-2404 


47 


127 


38 


9 


MACSJ2129.4- 


-0741 


82 


173 


51 


28 


MACSJ2214.9- 


-1359 


86 


209 


67 


19 


Total 




872 


2242 


627 


245 



about 4500 seconds. Bad-pixel masking, geometric distor- 
tion correction, cosmic ray rejection, image stacking and 
resampling were performed o n the flat-fielded image s us- 
ing the MultiDrizzle program (|Koekemoer et al.l l2002V For 
MACS J0018. 5+1626, observed in November 2010, addi- 
tional Charge Transfer Efficiency (CTE) corrections were 
applied t o the data using the p ixel-based CTE correc- 
tion code l|Anderson fc B edin 2010) before further reduction 
with the MultiDrizzle program. We applied the optimized 
resampling scales of 0.03", a Gaussian drizzle kernel, and 
pixfrac = 0.8 to avo id aliasing in the Point Spread Function 
l|Rhodes et al.ll2007h . Pseudo-color HST/ ACS images of all 
12 clusters are shown in Figure[T] Object detection and pho- 
tometry wer e performed on the redu ced ACS images using 
SExtractor dBertin fc Arnoutslll996r i. We removed the ob- 
jects with SExtractor FLAGS larger than 4 as well as point- 
like and spurious sources identified from their location in the 
distribution of peak surface brightness versus magnitude. 



3 ANALYSIS 

3.1 Cluster Membership Determination 

The most accurate way to determine cluster membership 
is via spectroscopic observations of galaxies in our target 
fields. Spectroscopically determined redshifts are available 
for ^ 1150 objects fro m C.-J. Ma et al. (2011, in prepara- 
tion) , iMaeTaD (|2008t ) and the NASA/IPAC Extragalactic 
Database. Galaxies with a radial velocity within 3<r (where 
a indicates the dispersion in the velocity distribution) of 
the systemic cluster velocity are taken to be cluster mem- 
bers (hereafter sample A), while those with a radial velocity 
above or below the 3<r threshold are selected as background 
and foreground galaxies, respectively. 

To increase the sample size, we alternatively select clus- 
ter members based on galaxy colour. For each cluster we 
identify the red sequence from a colour-magnitude diagram 
(F555W-F814W versus F814W; see Figure [2] for an exam- 
ple). The red sequence is formed by old elliptical galaxies 
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Figure 1. Three-color images of our 12 MACS clusters with F555W, (F555W+F814W)/2 and F814W used for the blue, green and red 
channel, respectively. Clusters are shown in order of ascending right ascension from left to right and top to bottom. In each cluster field, 
the red cross indicates the location of the peak of the X-ray surface brightness, while the yellow box marks the chosen BCG. The inner 
and outer white circles have a radius of 200 and 400 kpc at the cluster redshift, respectively. 
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Figure 2. Colour magnitude diagram of MACS J0717. 3+3745. 
The grey dots represent all galaxies in the field, and the black 
dots are red-sequence selected cluster members. The crosses mark 
galaxies with redshift information where green, blue, and red indi- 
cate spectroscopically confirmed members, foreground, and back- 
ground galaxies, respectively. 



whose spectra show similar Z?4ooo breaks resulting from pho - 
tospheric absorptions of heavy elements (|Bower et alJ fl992). 
Galaxies falling within ±2a of the red sequence (here a 
is the dispersion in the red sequence along the F555W- 
F814W axis) and featuring magnitudes of _F814W < 24 mag 
are classified as cluster members (hereafter sample B). The 
numbers of cluster members in each sample and for each 
cluster are listed in Table [1] 



3.2 Galaxy Orientation Measurement 

Galaxy orientations are measured with SExtractor from the 
F814W images. We use the windowed positional parameters 
from SExtractor with all derived quantities based on the 
weighted second-order brightness moments, 



Qi 



d 2 ew rg (\e\)e t e,i(e), 



(i) 



where W Tg (|#|) indicates a 2-dimensional Gaussian with 
a dispersion of r g , and 1(6) represents the brightness dis- 
tribution. The windowed positional parameters are derived 
adopting r g — FWHM/V8 In 2. The impact of the precise 
choice of r g is discussed in §4.3.2. We apply a Gaussian 
weighting function rather than constant weighting to avoid 
bias resulting from flux close to the noise level. Our analy- 
sis of galaxy orientations then uses the complex ellipticity 
e = e2 + ie2 and the object position angle, 8, as derived from 
the second-order brightness moments: 
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where 8 can be uniquely determined from the sign of 

Qxy (Fig. El). 



cluster 
centre 



Figure 3. Definition of 6, £, and <j>. 

3.3 Cluster Centre Determination 

The location of the cluster centre would ideally be defined 
as the position of the peak of the dark-matter distribution. 
Since, however, maps of the dark-matter distribution are 
not available for the majority of our clusters, we use two 
alternative estimators for the cluster centre: the location of 
the X-ray surface-brightness peak and the position of the 
BCG. We will further discuss the impact of the choice of 
cluster centre in some detail in §4.1.2. 

Seven out of the 12 clusters feature a single 
dominant BCG in the cluster centre; their aver- 
age offset between from the X-ray surface-brightness 
peak is ~5".7 (~35 kpc). The remaining five clusters 
(MACSJ0018. 5+1626, MACSJ0025.4-1222, MACSJ0454.1- 
0300, MACSJ0647.7+7015 and MACSJ0717.5+3745) con- 
tain multiple cD-type galaxies and show clear signs of 
substructure. For MACSJ0018. 5+1626, MACS J0454. 1-0300 
and MACSJ0647.7+7015, we choose the BCG to be the 
brightest galaxy within 200 kpc of the X-ray peak. For 
MACSJ0025.4-1222 and MACSJ0717.5+3745, we choose 
the brightes t galaxy in the mo st m assive sub s tructu re as 
identified bv lBradac et all I^POSt ) and lMa et all (|2009| ). The 
choice of BCGs are indicated in Figure [1] 



4 RESULTS 

4.1 Galaxy Orientation Distribution 

PK05 and F07 find that the major axis of member galaxies 
in nearby clusters tend to align with their radius vector. To 
examine whether such a radial alignment exists in our 12 
clusters at z > 0.5, we study the distribution of the values 
of <f>, defined as the angle between the galaxy major axis 
and the direction toward the cluster centre (see Fig. |3}. We 
stack the distributions in <f) obtained for all 12 cluster fields, 
using a metric scale to allow us to investigate the depen- 
dence of any alignment signal with distance from the cluster 
core. Note that BCGs are excluded from this analysis since 
their alignments reflect the surrounding large-scale struc- 
ture, as discussed in §Q] We use spectroscopically confirmed 
foreground galaxies as our control sample (hereafter sample 
C) since, unlike background galaxies that may be gravita- 
tionally lensed, they should not exhibit any preferred orien- 
tation. 

The distribution of galaxy orientations within the cen- 
tral 500 kpc (corresponding to the largest circle enclosed 
within the ACS images of all 12 target clusters) is shown 
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Figure 4. The distribution of galaxy orientation angles with re- 
spect to the radius vector; <f> = defines perfect radial align- 
ment. The black and red solid lines represent the distributions 
obtained when adopting the X-ray brightness peak or the position 
of BCG as the cluster centre, respectively. Error bars assume Pois- 
son statistics, and the dotted lines represent the average number 
of galaxies per bin. In the top-left corner of each panel we denote 
the respective subsample (A: spectroscopically confirmed mem- 
bers; B: colour selected members; C: spectroscopically confirmed 
foreground galaxies). 
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Figure 5. Same as Figure [4] but with the galaxy sample split 
according to distance from the cluster centre as indicated in the 
top left corner of each panel. 



in Figure [4] for samples A and B. We find no tendency for 
radial alignment for either sample; a Kolmogorov-Smirnov 
test (KS test) finds the observed distributions to be con- 
sistent with a uniform parent distribution at a probabil- 
ity of 31 (98)% and 62 (97)%, respectively, when the po- 
sition of the X-ray brightness peak (the BCG) is adopted as 
the cluster centre. As expected, the distribution of galaxy 
orientations for our control sample (C) is also consistent 
with random at a probability of 99 (99)%. Computing the 
average alignment angle ((f)) (where(0) < 45° indicates 
a preferential radial alignment), we find values of (if) of 
45°.45±1°.13 (44°.69±1°.15), 45°.21±0°.68 (44°.55±0°.69), 
and 44°.80±2°.26 (45°. 64 ± 2° .29) for sample A, B, and C. 

4-1.1 Dependence on Cluster Member Properties 

Although no alignment signal is detected for the largest sta- 
tistical sets of cluster members compiled for our sample, we 
need to bear in mind the possibility that radial alignments 
could be present for a subset of the cluster galaxies but are 
diluted to the level of insignificance when the entire galaxy 
sample is considered. Motived by the findings of F07, we thus 
examine if any radial alignment can be found as a function 
of galaxy colour or distance from the cluster centre. 

We first divide the samples A and B into two subsets, 
one including all galaxies within 200 kpc of the cluster cen- 
tre, the other comprising all galaxies within an annulus from 
200 to 400 kpc. A mild trend in favour of radial alignment 
is found in the central 200 kpc region for both spectroscopic 
and colour-selected cluster members when the position of 



the X-ray surface-brightness peak is adopted as the cluster 
centre (Figure [5} , yet a KS test still yields low probabili- 
ties of 77 and 78 % for the observed distributions deviating 
from uniformity. Similarly, the average alignment angle is 
measured to be (0) = 43°.43± 2°. 26 and 43°.47± 1°.24 for 
sample A and B. No trends at all are observed when the 
position of the BCG is adopted as the cluster centre. No de- 
viation from uniformity is detected in the annulus between 
200 and 400 kpc, regardless of how cluster membership is de- 
fined. The slight discrepancy between the results obtained 
in these two radial bins may be the result of a mild depen- 
dence of the strength of radial alignments on distance from 
the cluster centre. 

Using again the average alignment angle (<f>) as a mea- 
sure of deviations from a random distribution, we test for 
any dependence on galaxy colour by splitting our spectro- 
scopically confirmed sample of cluster members into two 
subsamples of red and blue galaxies, where "red" is defined 
as lying within 2a of the colour of the cluster red sequence 
(see H3.1[) . The results are shown in Fig. [5] (top panel) as a 
function of distance from the cluster centre. As noted be- 
fore, a mild trend of increasing radial alignments with de- 
creasing cluster-centric distances is observed, as evidenced 
by the fact that (4>) < 45° in the first three radial bins (bin 
size 100 kpc); however, only red galaxies contribute to the 
signal. Although statistically insignificant, these trends are 
qualitatively consistent with the findings of F07. 

In addition, we test for any dependence on galaxy lu- 
minosity. A similarly mild, but statistically more significant 
trend is found for the dependence of (<f>) on galaxy lumi- 
nosity (Figure [6] middle panel) in the sense that only the 
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Figure 6. Upper panel: (<j>) as a function of radius. Middle panel: 
(4>) as a function of absolute magnitude. Lower panel: absolute 
magnitude as a function of radius. The X-ray brightness peak is 
assumed as the cluster centre in this plot. In each panel, the whole 
spec sample of cluster members is in black, red and blue indicate 
red and blue galaxy members, respectively. 



most luminous red galaxies show a tendency for being radi- 
ally aligned. This trend is likely to be independent from the 
slight increase of {</)) toward the cluster core region, since no 
obvious correlation is found between a galaxy's luminosity 
and its distance from the cluster centre. 



4-1.2 Dependence on Cluster Relaxation State 

As a further test of whether subsets of our sample exhibit 
alignments that are diluted when the full dataset (Fig. 2} is 
considered, we examine the dependence of any signal on the 
relaxation status of the host cluster. Ideally, we would like 
to compare the alignment signal in virialized clusters with 
that found for morphologically disturbed, i.e., dynamically 
younger systems, and explore whether the strength of any 
alignments evolves as the host cluster approaches full relax- 
ation after a merger. To crudely quantify relaxation state we 
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Figure 7. As Figure|4]but showing the distribution of orientation 
angles for two subsets of sample A: the six clusters classified with 
Morphology code = 1&2 (top panel) and the six clusters classified 
with Morphology code = 3&4 (bottom panel). Both plots only use 
data from the central 200 kpc. 



use th e visual morphological classifications of lEbeling et al.l 
(|2007D which are based on X-ray morphology and the good- 
ness of the optical/X-ray alignment and result in a mor- 
phology code whose values range from 1 for fully relaxed to 
4 for extremely disturbed clusters. We here limit the anal- 
ysis to spectroscopically confirmed cluster members within 
the central 200 kpc, i.e. the data set for which a mild trend 
for radial alignments is found for the full cluster sample (but 
only when the X-ray peak is used to define the cluster centre; 
see Fig. [5] top panel). The results for clusters with morphol- 
ogy codes of 1 or 2 and 3 or 4 are shown in Figure [7] no 
significant difference between the two subsets is observed. 

Adopting an alternative, purely optical indicator of 
cluster relaxation state, we also examine the dependence of 
our findings on the presence of a single cD or multiple cD- 
type galaxies (see i]3.3|l . As shown in Fig. [8] (black lines), this 
trend originates primarily in the subsample of seven cluster 
featuring a single BCG, but again with very low significance 
(KS probability of 84%). The presence of just one or several 
BCG candidates makes no difference for the distribution of 
orientation angle when the position of the adopted BCG is 
taken to define the cluster centre (Fig. [8] red lines). 



4.2 Complex Ellipticity as a Measure of Galaxy 
Alignments 

Complementing our analysis based on orientation angles, we 
also search for galaxy alignments using the two-component 
complex ellipticity. Although the orientation angle and the 
two-component ellipticity convey similar information since 
they are mathematically dependent quantities, the two- 
component complex ellipticity has the advantage of allowing 
us to isolate a possible alignment signal while simultaneously 
testing for the presence of systematic effects in our measure- 
ment. Furthermore, the net alignment signal generated by a 
scalar gravitational cluster potential should be curl-free. Un- 
der the simplifying assumption that the cluster potential is 
isotropic, the curl-free and curl alignments can be decoupled 
by splitting the complex ellipticity e into a tangential/radial 
component and a cross-term component, i.e., 



Galaxy Alignments in Clusters at z > 0.5 7 




40 r 

30 h 




I 30 r 

20 [■ 








20 40 60 80 

cp (degree) 

Figure 8. As Figure[4]but showing the distribution of orientation 
angles for two subsets of sample A: the seven clusters featuring 
a single dominant BCG (top panel) and the five clusters with 
multiple BCG candidates (bottom panel). Both plots only use 
data from the central 200 kpc. 
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Figure 9. Variation of the two-component ellipticity with cluster- 
centric distance for Sample A. The solid and dashed lines repre- 
sent et and e x , respectively. Results based on the adoption of the 
location of the X-ray surface brightness peak (BCG) as the clus- 
ter centre are shown in black (red). To avoid cluttering, error bars 
are shown only for one set of results. 

e t = -S[|e| exp(-2i£)], (5) 
e. = -9f[|e|exp(-2ie)], (6) 

where £ is the angle between the galaxy radius vector 
and the image x axis ((/> = \8 — £[, see Fig. [3j. 

The variation of et and e x with cluster-centric radius 
is shown in Fig. [9] for the sample of spectroscopically con- 
firmed cluster members. We find a marginally significant 
trend toward negative e t in the central ~ 200 kpc region, 
consistent with the faint radial alignment signal detected 
in the orientation-angle distribution within the same radial 
range (Fig. [5] top panel) . No obvious trend is observed for 
the cross term e x , suggesting that systematic errors (other 
than the choice of cluster centre) are negligible. 



Figure 10. A comparison between the angular resolution and 
pixel sampling typical of SDSS (r' = 18.5 mag, left) and HST 
(F814W= 23 mag, right) images of galaxies. Note the difference 
in angular scale covered by either image. 

4.3 Possible Systematic Biases and Uncertainties 

4-3.1 Shape Measurements 

The high quality in terms of both angular resolution and 
depth of the HST/ ACS images available for our sample 
yields excellent pixel sampling and thus allows a much more 
robust shape determination than the groundbased data used 
in previous studies (see Fig. [TO] for an illustration). Es- 
pecially the spectroscopically confirmed cluster members 
(which tend to be brighter) are exquisitely resolved. In this 
section, we examine the robustness of our galaxy shape mea- 
surements under different weighting schemes (r g in §3.2) and 
detection thresholds. 

Throughout our analysis we use weights parametrized 
by r g = FWHM/V81n2 as provided by SExtractor. An al- 
ternative wei ghting scheme, giv e n by r a = FWHM, has been 
proposed by ISchrabback et al.l (|2007l ). In order to quantify 
the impact of the exact choice of r g on our results, we mea- 
sure the position angle of ~1000 galaxies in one of our cluster 
fields (MACS J2129.3-0741) for either value of r g . We find 
the average difference of (A9) — — 0°.13±1°.88 between the 
derived position angles to be negligibly small compared to 
the bin size used by us for the orientation-angle distribution 
(Figs. I4I5I7I8[) and conclude that our analysis is not sensitive 
to the details of the adopted weighting scheme. 

The impact of the chosen detection threshold can, in 
principle, be substantial but depends greatly on the weight- 
ing function. For the Gaussian weighting function applied by 
us in the shape measurement process, the resulting position 
angle varies by only about 1° when the detection threshold 
is raised from 5a to 20er. This effect too is thus negligible 
for the purposes of this study. 

4-3.2 PSF Anisotropics across the Field of View 

The anisotropy of the HST Point Spread Function (PSF) 
is known to be a potentially significant source of system- 
atic errors for studies involving galaxy shape measurements 
(e.g., weak-lensing analyses). The PSF varies across the ACS 
field of view and is also a function of time due to focus 
changes caused by "thermal breathing" of HST in orbit 
l|Rhodes et al.ll2007h . Since the cluster galaxies in our sain- 
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pie are well resolved we expect these PSF variations to have 
little effect on our study. We nonetheless test this expecta- 
tion by obtaining a quantitative estimate of the variations 
in position angle introduced by PSF anisotropics. 

We use MACS J2129.3-0741 as a test field since it con- 
tains ~80 stellar detections which allow us to sample the 
PSF across the field and fit position-dependent PSF mod- 
els. As the ACS/WFC images of MACS J2129.3-0741 were 
taken at six dithered positions with an integration time of 
~ 750 seconds each, the PSF model for this field will have to 
take into account both spatial and temporal PSF variations. 
We obtain a series of models for focus values ranging from 
— 10 t o 5 /im using the IDL Tiny Tim package ()Rhodes et al.l 
l2007f) based on the Tiny Tim HST PSF modeling software 
! Kristlll995l ). For each individual exposure we find the best 
telescope focus by comparing the model predictions with the 
shape of the stellar detections in the field, and then create 
and stack the PSF models for these six best-fit focus values 
according to the dither pattern used in our observations. 

Having thus obtained a PSF model for MACS J2129.3- 
07 41, we apply P SF corrections using the m ethod developed 
bv lKaiserl \ 19951 ) and lHoekstra et al.l (|l998h which measures 
the response in complex ellipticity to convolution with a 
small anisotropic kernel. The correction term in the com- 
plex ellipticity of individual galaxi es depends on their smear 
polarizability tensor, as defined in iKaiserl l| 19951 ). and on the 
second-order moments of the PSF model at each position. 
The resulting correction factor for the galaxy position angle 
is small though ((A6) = 1°.57) for the, in general, bright and 
well resolved galaxies in our sample. We conclude that PSF 
corrections have a negligible effect on the shape parameters 
measured for the galaxies in our cluster fields. 



5 DISCUSSION 

Evidence of radial galaxy alignments is sparse at best for 
the 12 very massive clusters at z > 0.5 investigated here. 
A marginal trend favouring radial alignment is found within 
the central 200 kpc region around the peak of the X-ray sur- 
face brightness, but becomes entirely insignificant once the 
uncertainties in the determination of the cluster centre are 
taken into account. Within a larger region encompassing the 
central 500 kpc, the 545 spectroscopically confirmed cluster 
members in our sample feature an average alignment angle 
of (<j>) = 45°38 ± 1°14, which is perfectly consistent with 
random. Our findings for clusters at z > 0.5 are thus in con- 
flict with the prominent radial galaxy alignments reported 
by PK05 and F07 for nearby clusters. Specifically, PK05 
find {(j>) = 42°79±0°55 for 2200 spectroscopically confirmed 
cluster members within the central 2 Mpc of the centres of 
clusters at z ~ 0.1. We note that this discrepancy can not 
be attributed to statistical uncertainties: the signal reported 
by PK05 would still be significant at the 2<r confidence level 
if their sample size were reduced by a factor of 4 to 550 ob- 
jects. The complete lack of evidence of a net radial galaxy 
alignment for our z > 0.5 sample within a much smaller 
distance from the cluster centres is particularly puzzling if 
radial alignments are indeed strongest near the cluster cores, 
as suggested by both our results and F07. In this section we 
discuss differences between our work and previous studies 
that might explain the discrepant results. 



We first note that the mass range probed by our clus- 
ter sample is considerably different from that used by F07. 
MACS, by design, selects the most X-ray luminous clusters, 
and the ones selected for this work feature total masses in 
excess of 5 x 10 14 M^ jBradac et alj|200gl ; ISmith et al.ll2009l : 
iMantz et all l201Ct E imousin et al.l 120111 ); by contrast, the 
sample compiled by F07 using a group finder algorithm cov- 
ers a much wider mass range of 5 x 10 12 Af© < M V i r i a i < 
5 x 1O 14 M0. The clusters in the sample of PK05, however, 
are X-ray selected just like ours and tend to be more massive 
than those studied by F07, and yet these two studies find 
consistent radial alignment at low redshift. We conclude that 
the lack of alignments observed by us at z > 0.5 is unlikely 
to be caused by a difference in cluster masses. 

Another possible systematic difference between ours 
and previous work is the completeness of the spectroscopic 
follow-up observations of cluster galaxies. Spectroscopic cov- 
erage tends to be inhomogeneous across clusters due to con- 
flicts between slits or fibers specifically in the dense cluster 
cores. Since the strength of any net radial alignment of clus- 
ter galaxies may depend on their distance from the cluster 
centre, spectroscopic completeness as a function of radius 
needs to be taken into account when comparing results ob- 
tained for different datasets. Extensive spectroscopic follow- 
up observations of our 12 MACS clusters result in an overall 
completeness of about 80% down to F814W= 21.5 mag. As 
shown in Fig. [11] the spectroscopic completeness above this 
limiting magnitude varies by only ~ 15% from the central 
200 kpc to the edge of the ACS field of view. Correcting for 
relative incompleteness of the spectroscopic coverage does 
not change our finding that no significant radial alignment 
is detected within the central 500 kpc. Although a compar- 
ison with PK05 is not possible since no information of this 
kind is provided for their sample, we cannot imagine a selec- 
tion bias in the spectroscopic survey underlying the PK05 
dataset that would explain the stark discrepancy between 
their results at z ~ 0.1 and ours at z > 0.5. 

Differences between our results and those in the liter- 
ature could also be caused by biases affecting the measure- 
ment of galaxy shapes and orientations. We examine pos- 
sible systematics in our analysis of HST/ACS data in §4.3, 
and find that both instrumental and measurement system- 
atics are negligible. Instrumental systematics should also be 
negligible for the datasets used by PK05 and F07 since the 
SDSS is conducted in drift scans, paying multiple visits to 
each location. Systematic errors, however, can be expected 
to be relevant for shape measurements based on ground- 
based SDSS data, owing to poor seeing and insufficient pixel 
sampling. 

Other potential biases inherent to alignment studie s 
based on SDSS data are investigated bv lHao et al.l (|20 1 lh - 
Using SDSS DR7 data, they find that the radial alignment 
signal correlates with the apparent magnitude of BCGs but 
not with their absolute magnitude, causing them to conclude 
that such alignments are created by systematic errors in the 
galaxy orientation measurements caused by contamination 
from diffuse light of BCGs. While such a bias appears plau- 
sible in principle, it remains unclear how diffuse light from 
the BCG can affect the measurement of galaxy orientations 
out to the large radii over which highly significant radial 
alignment is reported by PK05 and F07 (2 Mpc and at least 
0.5-Rviriai, respectively). Furthermore, since the control sam- 
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Figure 11. Spectroscopic completeness of cluster members with 
"2F814W < 21.5 as a function of clustercentric radius for our 
sample of 12 MACS clusters at z > 0.5. 

pie used by PK05 consists of field galaxies in cluster fields, at 
least some radial alignment trend should be present in their 
control sample too, if the measurement is indeed biased by 
diffuse light from BCGs. To summarize, the systematics af- 
fecting the results obtained from groundbased observations 
are still controversial and need further investigation. 

Finally, a possible physical explanation for the differ- 
ences in alignment strength between our z > 0.5 sample 
and samples comprising nearby clusters is that radial align- 
ment evolves dramatically with redshift. If radial alignment 
is the result of interactions of cluster galaxies with local tidal 
fields, cluster members newly accreted from the field pop- 
ulation may not have been orbiting the cluster centre long 
enough for tidal torquing to have a measurable effect. This 
effect can be expected to be stronger in younger clusters, 
i.e., at higher redshift. Indeed, clusters at higher redshift are 
known to be less relaxed than local ones (M arm & Ebeling 
Although we find no correlation between cluster mor- 
phology and radial galaxy alignments within our sample 
(see i)4.1.2[) . perturbations from minor or major mergers ex- 
pected to be more common and frequent at larger lookback 
times may well prevent a radial alignment signal from reach- 
ing a significant level in distant clusters. 



location of the BCG is chosen as the cluster centre, the sig- 
nificance of any deviation from uniformity in the alignment 
angles is even lower. We find no difference regarding galaxy 
orientations between two subsets of our cluster sample de- 
signed to separate the more relaxed from the more disturbed 
systems. 

In order to test for systematic errors in our analysis, we 
examine the distribution of both components of the com- 
plex ellipticity; as expected for curl-free alignment signals 
generated by a scalar gravitational cluster potential, the 
cross term vanishes at all clustercentric radii. We also verify 
the negligible impact of systematic uncertainties and biases 
caused by the variability of the HST point-spread function 
(both temporal and across the ACS field of view), and by our 
choice of weighting function and detection threshold when 
measuring galaxy ellipticities and orientations using SEx- 
tractor. 

The absence of significant radial galaxy alignments in 
our 2 > 0.5 clusters stands in stark contrast to the promi- 
nent radial alignment reported for nearby clusters. Specifi- 
cally, we find no significant radial alignment of cluster galax- 
ies within the central 0.5 Mpc for clusters at z > 0.5, whereas 
a strong alignment signal is detected out to 2 Mpc for clus- 
ters at 2 ~ 0.1. We examine several systematic effects that 
could cause this discrepancy and rule out differences in the 
selection of clusters or cluster members. The most plausible 
explanations are (1) that the radial alignment evolves dra- 
matically with cluster redshift, in the sense that the signal 
is weaker in dynamically younger systems, or (2) that the 
signal observed in nearby clusters is at least partly spurious 
and the result of measurement biases in the low-resolution 
SDSS data used by the respective studies. Both of these ex- 
planations are testable and will be examined in a future, 
extended study. 
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6 SUMMARY AND CONCLUSIONS 

Using high-quality HST/ ACS imaging data, we have per- 
formed a search for galaxy alignments within the galaxy 
population of 12 very X-ray luminous clusters at z > 0.5 that 
constitute a statistically complete subset of the MACS clus- 
ter sample. We find no evidence of a net alignment within 
500 kpc of the cluster centres ({<(>) = 45°38 ± 1°14, based on 
545 spectroscopically confirmed cluster members). Within a 
radius of 200 kpc of the cluster centre (defined as the loca- 
tion of the peak of the X-ray surface brightness), the distri- 
bution of galaxy orientations shows a mild, but statistically 
insignificant trend favouring radial alignments (according to 
a KS test, the probability of the observed distribution be- 
ing different from a uniformly distributed parent sample is 
less than 80%). This slight alignment trend is found to orig- 
inate in red galaxies near the cluster centre. If, instead, the 
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